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ABSTRACT 

O . . 

^vq ' We give a rigorous treatment on the foundation of the first order asymptotic theory of quantum es- 

^->' timation, with tractable and reasonable regularity conditions. Different from past works, we do not use 

^ ■ . . . 

Fisher information nor MLE, and an optimal estimator is constructed based on locally unbiased estimators. 



^ 






o 
o 



Also, we treat state estimation by local operations and classical communications (LOCC), and estimation of 



^^ ' quantum operations. 



1. INTRODUCTION 

The purpose of this paper is to give a rigorous foundation of the first order asymptotic theory of quantum 
estimation, which has been established in these years. In addition to most basic setting, we also treat state 
estimation by local operations and classical communications (LOCC, in short) and estimation of quantum 
,,_^ ^ operations. 

0^ , This research field was initiated by Nagaoka(1987), Nagaoka(1989), followed by Hayashi and Mat- 

[■^~. . sumoto (1998), Gill and Massar (2002). (Many of important papers in the field are included in Hayashi (2005).) 



^^ ' optimal efficiency of consistent estimators, they had reduced the optimization of consistent estimators to 



optimization of Fisher information, or equivalently, of locally unbiased estimators. These works had laid 
foundation on which number of works, mostly computation of asymptotically optimal estimators and their 
/\^ • costs, are based. In closer look, however, they either miss the detail of the proof, or assume intractable 

Cu ' regularity conditions. 

One reason for such incompleteness is that the focus of these works were consequences of the foundations, 
rather than their rigorous proof. Also, the following technical difficulties seems to be a part of reasons. In 
quantum statistics, the probability distribution of the data depends on the choice of measurement. Therefore, 
for classical estimation theory to be applicable, a set of regularity conditions should hold for all the probability 
distributions resulting from arbitrary measurement of interest. In Hayashi and Matsumoto (1998), they use 
this sort of statement as their regularity condition. As a result, their regularity conditions are quite difficult 



to check for given quantum statistical models. 

The purpose of the paper is to provide rigorous proof assuming tractable regularity conditions, including 
the case of infinite dimensional Hilbert space. In addition to the most basic settings, we also treat state 
estimation by semi-classical measurement and by local operations and classical communications (LOCC, in 
short). Also, estimation of quantum operations is studied. 

Different from previous works, we avoided use of Fisher information, and composed an asymptotically effi- 
cient estimator from an optimal locally unbiased estimator, because of the following reasons. First, quantum 
asymptotic Crammer- Rao bound is not a simple function of any quantum analogue of Fisher information. It 
equals Holevo bound, which is defined in terms of operator version of asymptotically unbiasedness conditions 
(Hayashi and Matsumoto (2004), Matsumoto (1999), Guta and Jencova (2006)). The second motivation is to 
simplify the regularity conditions, by avoiding technical difficulties stated above. 

One of major difference between quantum mechanics and classical mechanics is behavior of composite 
systems. In quantum mechanics, the state of the system and the measurement in composite systems may 
not be in convex combinations of those without correlations between subsystems. In such cases, we often 
observe non-trivial quantum effects, which can never be reproduced by classical mechanical random variables, 
such as violation of Bell's inequality. Therefore, it is of interest to compare measurement with non-trivial 
correlations and the one without it in their efficiency of state estimation. 

For that purpose, we study semi-classical measurements and LOCC (, short for local operations and 
classical communications,) measurements. In the former, we are not allowed to use measurement collectively 
acts on given n independent samples. In the latter, each sample is a state in a composite system (A and B, 
say), and we are not allowed to use the measurement quantumly correlating over A-B split. 

The last topic is estimation of a quantum operations. It had been observed that for some cases (e.g., 
unitary operations, or noiseless operations), the mean square error of optimal estimators scales as O (l/n^j 
(Heisenberg rate), which is significantly smaller than 0(l/n), and there had been suggestion of efficient 
measurement scheme utilizing this effect. Recently, however, several authors ( Fujiwara(2005), Zhengfeng 
Ji, et. al. (2006), etc) had pointed out that O (l/n^)-scaling is not observed in some class of operations 
(typically, they corresponds to noisy operations). We show that O (l/n^)-scaling is rather exceptional, and 
not observed so long as the model lies in interior of the totality of quantum operations. 

2. QUANTUM ESTIMATION THEORY 

2.1. QUANTUM STATE AND MEASUREMENT 

In quantum mechanics, the probability distribution of data z G M' is a function of the state p of the system 
of interest, and the measurement M which is applied to the system. The probability that uj lies in a Borel set 
A, the corresponding random variable, and the post-measurement state is denoted by P^ (A), fi, and p^, 
respectively. (Throughout the paper, the random variable is denoted by capital letters, and the elements of 



its range is denoted by its dccapitalization.) 

p and M are represented by linear operators defined in a separable Hilbert space ("H, say). The inner 
product of (p and tj: is denoted by (iy9, ■0). We assign to the composite of the system "Hi and 'H2 the tensor 
product "Hi ®'H2i which is the linear span of {ei^i ® &2,j} ( {ei,i} and {£2,1} be a complete orthonormal basis 

(CONS) of -Hi and -H2 respectively). 

1 /? 
The notation \A\ means \A\ :— {AA^^ , and ||A||j^ :— trjAj is a quantum version of total variation. 

The totality of trace class operators, or operators with ||A||-|^ < 00, is denoted by tc{'H). Also, ||j4|| := 

sup||„|j^i ||A(^|| and B{'H) denotes the totality of hounded operators, or operators with ||A|| < 00. (the 

standard norm in M™ and in % is also denoted by ||-||.) We introduce an order in the space of matrices 

hy A > (>) B <^ {(fi, Aip) > (>) {if, Bif), Vip. An operator A is said to be positive, if A > 0. A mapping 

A of Tc{H) to Tc{H') is called completely positive, if A ® I : B{T-L®JC) -^ B{T-L'®JC) is positive, i.e., 

A > ^ A ® 1(A) > 0. A is said to be trace preserving if tr A = trA(A) (VA). Also we define 

l|A|L6:-sup;f^l|;f|l^^i||A®I(A)||i. 

A state of the system is represented by a density operator, or an operator p with p > 0, p ~ p*, and 
tr/3 = 1. A measurement M is represented by an instrument, or a cradditive map M : A — > M [A] of the 
collection *8 of Borel subsets in M™ into a completely positive linear transform M [A] in tc {%) with M [K'] 's 
being trace-preserving. Here, a-additivity is in the sense of strong operator topology in B {tc{'H)). Using 
p and M [A], P^ (A) and pA is given by tr M [A] p and -pwr^^ [^] (p), respectively. An operation which 
does not extract information is described by a completely positive and trace-preserving (CPTP) linear map 
A from Tc(-H) to Tc{n'). 

When we are interested only in Pj^ (A), we use a positive operator valued measure [POVM, in short), or 
a CT-additive map M : A — > M (A) of S to positive Hermitian operators with M (R') — 1. Here, cr-additivity 
is in the sense of weak operator topology in B{H). The POVM M corresponding to the measurement M 
satisfy P^ (A) = tr M [A] p = tr pM (A). Throughout the paper, POVM of a measurement is denoted by 
the same character as the measurement but in the standard font. 

The support supp (M) of the instrument M over *B (R') is the smallest set with M [supp (M)] = M [W] . 
The support of a POVM and a measure over *B (R') are defined analogously. 

In this paper, we need integral of the function taking values in tc {%) and B (tc {%)), which is a Banach 
space with the norm ||-||j^ and ||-||^t,, respectively. A Banach space valued function / is called strongly 
measurable iff Ve > 3/' ||/ (x) — /' (x)|| < £ holds almost everywhere. / is called weakly measurable iff 
{y* , f (x)) is measurable for any element y* of the dual space. Since tc{'H) and B{tc{'H)) are separable, 
these two concepts are equivalent in our case due to Theorem 1.1.4 of Schwabik and Guoju (2005). 

Pettis integral of weakly measurable function / is defined by the relation J {y* , f (x)) dx — {y* , J f (x) dxj, 
Vy. Bochner integral of a simple function "^iCiXAi is defined as ^jCi/i(Ai). For a strongly measurable 



function /, it is defined as lim„_>.oo / /n (x) dx (convergent in norm), where {/n}„ is a sequence of simple 
functions with lim„_i.oo \\fn (x) — /(2:)|| = ahxiost everywhere. Bochner integral exists iff J ||/|| da; < oo 
(Theorem 1.4.3 of Schwabik and Guoju (2005)). Fubini's theorem holds for Pettis integral and Bochner 
integral. 

2.2 ASYMPTOTIC THEORY OF QUANTUM STATE ESTIMATION 

Suppose that we are given n independently and identically prepared samples, i.e., the system "H® ■ ■ ■ (gj'H :~ 

n 

■H*^" in the state pe (d ■ ■ ■ ® pe ='■ Pg , where pg is drawn from a quantum statistical model A4 := {pg ; 9 G 0}, 

n 

with 0's being an open convex region in M™. 

Our purpose is to estimate the true value of 9, based on a measurement M" acting in "H**". Based on the 
measurement result Wn e M'" , we compute the estimate r„ of 6*. The pair £„ := { M", T„} (or sometimes the 
sequence {fn}^i also) is called an estimator. T„ is a measurable function of M'" to 6„ C M™. The follow- 
ing notations are used: E^^" [/ (w„)] := / / (w„) trp^A/" (dt^„), (MSEe [En]),^^ := E^^" (T^ - 9') {T^ - 9^), 
(Ye [£n]),^j ■■= K" {Tn - K^^ [Tn]) (^^ " K'^ [T[\) ■ Below, Tr denotes the trace over R™, and 9, := -^. 
Gg is a symmetric positive real matrix, and 9 — > Gg is continuously differentiable, TrGe < &i, and 
iTrGe-TrGs'l < hi\\9-9'\\. We also define {Bg^[£n])] ■= SjE^^^" [r^]|g^g^. Our interest is the first 
order asymptotic term of the weighted mean square error lim„_i.oo riTrGgMSEe [£„], minimized over asym.p- 
totically unbiased estimator , or {£n}n=i with the following condition: 

lim E^ [Tn] - 9, hm {Bg [£^])] = 6), V0 E 9. (1) 

In considering ([T|), Eg^ [r„] has to be differentiable, which is made sure by Lemma[2J Use of MSE 

may be justified based on the existence of the asymptotic normal efficient estimator, which is composed in 

Subsection 3.3. 

Our purpose is to replace this condition by the following tractable condition without changing the optimal 

lowerbound to the asymptotic cost: £gg_n — {Mg , Te^^n} is said to be locally unbiased at 9q if 

<'Mreo,«] = ^0, {Bg„\£g,,r^)]^6]. (2) 

Note that the condition ([2]) is closed at the point ^o- In the following sections, we prove that minimization 
of lim„_>.oo ?^Tr GgMSEg [£„] over all the asymptotically unbiased estimators can be reduce to minimization 
over the locally unbiased estimators under some proper regularity conditions. 



3. THE BASIC SETTING 

3.1. REGULARITY CONDITIONS AND ASYMPTOTIC CRAMER-RAO BOUND 

Regularity conditions on quantum statistical models and estimators are listed in Table[l] in which convergence 
is with respect to ||-||^. Oi,e,n is as defined in LemmaO and 0,- g „ := 0i,e,n®pf"~"^+Pe®0i,e,n'8'pf "^^H • 



(M.l) dipe and didjpe exist and are locally uniformly continuous. ||9jp0(,||,||9i9jpe||j^ < oi < oo. 
(M.2) 3Le,i: Hermitian and dipe = ^ {Le.iPe + PeLe.i), and ti pg {Lg,i) < oo, V6' e O. 
(M.3) There is an estimator £„ = |m", fnj in "H®" , such that 
(M.3.1) dl]) and (E) are satisfied. 



T„~e 



< ^, yeee.3 d. 



(M.3.2) E 

(M.3. 3) M" is n times repetition of a measurement M in %, producing the data x^, G 



(E) 3a4,„, V0 e e, / ||r„ K) - 9\\ tvOl'^M" (da;„) < naia4,„, / ||T„ K) - ^f trpf"M" (dw„) < nal „. 



(E') T„ takes values in Ot„, with supg g,^Q \\9 — 9'\\ < a^^n < co. 



Table f : Regularity conditions on quantum statistical models (M.f-4) and estimators (E), (E') 

Among the conditions on models, only (M.f) is needed to prove the lowerbound. Unless otherwise 
mentioned, (M.l) are assumed throughout the paper. (M.2-3) are necessarily to prove the achievability of 
the lowerbound. (M.2) is equivalent to \diti- pgX\ < cltr pgX^\ for any bounded Hermitian. 

If dim H < cx), an example of estimator £„ = -^ M". Tn > with (M.3. 1-3) is constructed as follows. Let I := 

'J^ 7 1 

(dim■H)^ and define e„ := fo, • • • ,0, 1, 0, • • • , o") G R'. Let supp fM^) he {e^}l^^, and let |M({e„})| 
be linearly independent. Denoting the K-th measurement result by wi_«;, we can estimate tvpgM {{ey}) by 
the relative frequency of observing e^, which is i7-th component w^ of wi :— — X]"=i ^i,k,- Let p be a solution 
to the system of linear equations trpM ({e^^}) ~uj^ {v ~ 1, ■ ■ ■ , I), and T„ is defined by p^^ = H (p), where 
n is a properly defined projection. Also, if {pe}g^Q is a smooth submodel of quantum Gaussian model {cr,,}, 
we can compose £„ based on the estimator f)„ of 77 by p^^ = H (o'^„), with proerly defined projection H. 

Both of them has the following property. {pg}g^Q is a somooth submaniforld of a larger quantum state 
model {o'jj}, where ry has consistent estimator in the form of f/n = ~Sk=i'^1'«' where wi^^ is the data 
obtained by application of M on the K-th sample. Suppose that 77 = (0, C)j ^nd pg — (TgQ(gy Moreover, we 
suppose that C, (0) is uniformly continuous in 6. Then, T„ :— {fj\, ■ ■ ■ ,i]™) satisfies the requirements. 

As for the estimators, besides ([Ij, we suppose f„= {M", r„} satisfies (E) in Table[2]for all n. (E') is used 
to characterize lowerbound to the asymptotic cost. Observe that (E')=>(E). 

We define the asymptotic quantum Cramer-Rao type bound Cg {Gg,A4) as 
lim inf {nTr GeMSEe [£„] ; M" in H^", ([1]), (E)}. In the succeeding subsections, the following theorem will 

n— ^00 

be proved. In the remaining of this subsection, some technical lemmas will be shown. 

Theorem 1 Suppose (M.1-3) hold. Then, 

C^{Gg,M)^ lim inf{nTrGeV(,[£e„] ; M" m'H'^", 0j, (E') |, (3) 

= lim mi{nTTGgYg[£e,n] ; M" m'H'^", 0), (E) }. (4) 



Lemma 2 (E) and (M.l) imply the existence of djE'^" [T^] and djE'^" [T^] = /T,* (lj„) tr^^peM" (dw„). 

Proof. Due to Lenima[21 this Lemma is equivalent to Proposition VL2.2 of Holevo (1982). ■ 

Lemma 3 (M.l) implies that 3ai3a2Vi, V9, 9^ G Q and \9^ ~^o| < '^2,- 9^ = 9q (j =^ i), 30i,e such that 
|9iP0ol < Oi,e , trOj,e < ai < oo. 

Proof. Since dipe^ = dipe-a^e, + SxLe-a2e, ^fp^'^^^ '^^■fi '■= \d%Pe^a2e,\ + Sx^eTa^^, I'^fP^;! dx, if exists in the 
sense of Bochner, satisfies requirement. This is true since 9j^p6» L is continuous in 9 (hence, measurable and 
integrable over the finite interval). ■ 

Lemma 4 (E'), combined with (M.l), implies 

a*^a*''Er [Tr,]^ jT^{iOn)trdfdl'pf^Ar{doJn) {t,,tk e {0,1}), 



TrGeVe [En] - TrGe-Ve- [f„]| < [na^ + 1) 6i (a4^„)' \\9 - 9'\\ , 
E^" [T„] - E^," [r„]|| < m2na4,„ai \\9 - 9'\\ , 
a,-E^ [r„] - 9,E^; [r„]|| < m^n^a^^naf \\9 - 9'\\ , 
lim {Bog [£e,n])] = Sj, where {■^eo.nleose satisfies f^i- 



(5) 

(6) 
(7) 
(8) 
(9) 



Proof. (E') implies | / T^ (w„) tr tM" (dw„)| < |/ JT^ (w„)| trrM" (da;„)| < ||r||;^a4,„. Therefore, the map 
r — )■ / T",* (wn) tr tM" (dw„) is a continuous linear functional, and is exchangeable with lim. Therefore, the 
first two identities follow. To show (O, apply the mean value theorem to the function 9 -^ E^ [T„]. Due 
to (HD, we obtain |E^" [T^] - E^," [T,^] | < ^^^li \JT^ {^n)tvdjpf^W (da;„)| \9^ - 6I'J|. Therefore, due to 
(M.l) and Lemma[3l we have ([7]). ([U is shown similarly. To show ([9|), observe 

{Bea[£0,n])]- 5] = (Se^ [fe.n])* - (-8^0 [fe(,,„])^ 



< 



{Be,[8g,,])]-{Be[8g,,])] + (B^ [fe,™])^- - (Be„ [f^o,™]) ■ = {Bg^lEg,,])] ~ {Bg[Eg^n]) 



Due to dl]), we have Q. 



3.2 LOWERBOUND AND ^^® 

First we prove that the RHS of (|3]) is a lowerbound to Cg {Gg,Ai). Define locally unbiased estimator 
£g^„ = { M", Tg^„ } by r„ = Be [£„] {Te,n -9)+ E^^" [T„]. Obviously, 

nTr GeMSEe [£„] > nTr GgYe [En] = nTr GgBg [£„] Ve [f e,«] S^ [£„]^, 
and letting rt — >■ oo, we have our assertion due to ([T]). 

Below, we prove (l3]) = (l4]). Since (E') implies (E), it suffices to show ([3])<(|4]). Suppose £g^n satisfies (E) 
and Q. Let S"^,, := Tg.n in |jr0,„ - 6*11 < L-case and S^^ := 9 otherwise. Let J'|;„ := {M^,rg^„|, and 



6 



let Eg.^ = (Me,T'/'„| be a locally unbiased estimator with T^,^ = Bg J^^^ (5^„ 



-E 



^a 



Obviously, £q^^ satisfies (E'). Also, due to LemmalHl Lemma[31 and (E), we have 



[Be [P^ ) \ - <5} = dj (tI,^ - 0^) P;« (d..„) 



Te,,^-e\\>L 



Te,„-6I||>L 
< / \\Tg,„-0\\tT0\'^^M^{diOn)^0 (L^C»). 

To„-e||>L ' ' 



n^-eAtrd.pf^^M^idiOn) 



s(l) 



Therefore, Ve > 03L, 
TTGeVe\£e,n]>TrGeYe 



n 



TiGeBg 



n 



Ve 



Bg 



•'a 



>TrGeVe 



Taking infimum of the both ends, wc have ((3])<(l4|). 



3.3 ACHIEVABILITY 

Based on {£e.ni}eeQ — {M^^ , ^e,ni jeee such that ([2]) and (E') with n ~ ni are satisfied, we construct a good 
estimator £^^1 with 2 steps in the following. Given p^"" , invest p^"" to obtain the data wi := (cjii,--- ,cji,„q), 
where LUi^i G M'. Based on the data, we compute the estimator 6*0 = T„(, (wi). Now, we divide pf"'~"'° into 
the ensembles each with rti copies. The munber of ensemble, "~"° , is denoted by n2- Here, no and n2 are 
chosen so that hq — n^ is satisfied. We apply M^^ to each ensemble pf''^ obtain the data W2,i, • • • , 
W2,n2(^ M'"i) and compute 



H "2 

r"' := V^T'so^nj (CJ2,; 



(10) 



K = l 



The measurement defined above is denoted by M*" 



Lemma 5 Suppose that (M.1,3) hold. Suppose also that the family {fe.niJsGe satisfies (0) and (E') with 
n — ni, V0 G O. Then £^^ constructed above satisfies lim nTrGeMSEg [£"^] < ni lim TrGeVg [i^eo.ni]- 

122— >oo 6*0— >9 



Proof. Applying mean value theorem to the function 6* — >■ E^ 



Co 



r, 



00, "1 



we have 



\ / u — Oq 

(11) 

where 7^^* is the reminder term. With the help of ([5]) and (M.l), 



')'e,Vn 






j,fc=i 



<nfm"afa4,„J|eo-e|r, (12) 



where 9' lies between 9o and 6*. Since MSE is the sum of the variance and square of the bias, we have 



TrGeMSEe 



^n' I Trig - 0a 



= TrGeVe 



CMT„o =0o\ + Y. (G«),,7,"V:7,"V' 



jj=l 



< — TrGeVe [fe„,„J + ntm^ (a?a4,nj' TrGe ||0o - ^11"* 

7l2 



Taking average over T„q of the left most and the right most end, 



hm TiTrGeMSEg [f,;^i] < hm 

712— ^oo n2— >-oo 



-E,^"°TrGeVe 

"-2 



F - 



ifrTi"* {a\a4^^ni) TrGgE, 



M"o 



Tno - 



-^-^ sup TrGeYe [£eo,ni] + nnjm'^ (alai^m) TrGgEj^ ° 



Tno - 9 



< hm sup mTrGgVe [i^eo,"!] 
(i)"2^°°eo:||eo-e||<e ^""-o Sos 

< hm sup TiiTrGeVe [fso.ni] H 77T(a4,ni) TrGe + hm (712^1 + no) "•i^"' (0104 „ J ^r-bTrGe 

rii^ n2— >-00 /] .11/) flll^.^ (r4/Ki^ n2— >00 'xi ' 

U^J t'OnlPQ — y|| <£ C TI2 TI2 



= sup niTrGgVg [£9o,n^] ■ 

eQ:\\eo-e\\<e 



Here (i) is due to P^"° < T„g — 9 > s> < -^-^ which fohows from (M.S. 2) and Chebyshev's inequahty, and 
(ii) is due to (M.S. 2). Since e > is arbitrary, the lemma holds. ■ 

Lemma 6 Suppose that (M.1-3) hold. Then {£n^}^^i satisfies (E) and (Qp. 



Proof. Observe \\T^^ - e\\ < 
satisfies (E), also. 
Observe 



Tno - S 



a^^n holds. Since T„j, satisfies (E) due to (M.3.1), {f"^}^i 



ES 






<E 









f m 



(ii) 



M"0 



<E, 

(i) 
2 



M"o 



So,T„a 



,2^2, 



< nfm^aia4,„iWE^^^ 

(iii) 



M"0 



Tno ~ ^0 



(iv) 



Tna — ^0 

Here, (i), (ii), (hi), and (iv) is due to pT|) . p^ . concavity of y^, and (M.S. 2), respectively. Therefore, 
-gM"i" [t^u] —^ 6»j^. (Sg^ [fn])* -)■ (5* is proved as follows. In Subsection 4.S right after the statement of 
Lemma fTTl we will prove 



ddE^ "E« 



T;i 



8=00 



= d,, E^ "E 






Ti 



Er"a:E/- 



0=60 



Ti 



e=eo 



(IS) 



Defining L",- := La ,;®1^ 



i+l®L0,,®l®"-2^- ■ ■+l«5"-i®Le,, , we have 9,pf " = i (i^,,pf " + pf "^e,.) , 



trpf" L^\J =ntrpe(Xe,z) , and 



a»trpf"A = tra.pf"A ==mrpf"AL^„ VA: bounded Hermitian, (14) 

where the first identity is due to the continuity of linear functional X — > tr XA (e.g., TheoremH.7.2 of 
Holevo(1982)). ([Ti|) . in combination with Schwartz's inequality, leads to |9itrpf"X| < nti pg (Lg^i)'^ tr pf'^ti X"^. 



Observe 

we have d. 



7fl 



E/"° 



Ti 



0^ 



< 



T - 0^ 

-i-na ^0 



04. „j. Hence, due to (fT2|) . (M.S.I), and Lemma[2l 



M"o„ni J 



E,-' 7, 



»0,-<,io 



/ 7^'i tra,pf^"''M"o (dcJi). Therefore, due to Theorem VI.2.1 of Holevo (1982) 



and (|14p . the first term of (|13p is evaluated as follows (they are used to show (i) below). 



9,;(E«^""E, 






5,, E 



nM"o 



n3 I ^"iJ 



ft 



W '0„ T 



< 



:i) \/trprK.)'<"te„o)' I) \l-otrpeALe..f- 



2 nim^aia4,„j 132,60 



where (ii) is due to (|12p and (M.S. 2). Therefore, the first term vanishes as no — >■ oo. Due to (j9]) of 



LemmaSl the second term converges to ft Eg " 






e=eo 



(5*, and (P) is proved. 



Lemma 7 lim inf {TrGeVe [feo,n] ; &, (E')} = inf {TrGeVe [8e,n] ; &, (E')} 

©o— >6 

Proof. Suppose the LHS is larger than the RHS (, denoted by A in the proof) by 2c > 0. Then one can 
find a sequence {Ok} such that lim inf {TrGgVe [£9^,11] ', ©J (E')} = ^ + 2c. We prove this cannot occur. 

Obviously, among those satisfying ([2]), (E'), one can find {£e,n}g^Q such that TrG gVg [£e.n] < ^ + 
c. Define f ^^ „ := {m^, T^^_„} by r^^_„ := Be, [^e,™]"' (t^,,, - E^J [re,„]) + 9^. It is easy to verify 



Ve 



£, 



-Be, [Se.n] ^ TrGeVe [i?0,n] (-60, [Se.n] ^) and that f g^ „ satisfies ([2]) and (E'). (Here note that 



a4.„ has to be replaced by the other constant.) Therefore, due to © of LemmalU 



lim TrG eVg 



Sfc.n 



= TrGeVe [£e.n] <A + c<A + 2c= lim inf {TrGeVe [£e,A ; ©, (E')}. This is 



contradiction. ■ 

Due to Lemmas[S][71 we have '<' of ([3]) of Theorem[TJ 



3.4 On asymptotic normality of the estimator ()10p 

The estimator (jlOp is asymptotically normal. We prove the assertion in ?n = 1-case, supposing that 

infeogn™ Ye ['?9o,"i] is not 0. 



Pr"" {V^V, [f,,„,]-^ (n^^ -e)<y]-^ {y) 



E^°R 



<E 



M„o 



R ""« ^AIV, 



T„(,,ni 



(rpm 



'e.T„ 



■E 



M„o 



IVI - -. j^ -. IVI -. I 

p, ""» {V^Ve [£(,,„,]-^ (t;'^ - 0) < yj - p, """ V^v, 



° {V^V, [£e,„J^^ {t:^ - 0) < y} - $ (y) 
< y 1^ - * (y) 

^(i:'^-^-7;x)<y 






Due to Berry-Esseen bound (Chapter 11 of DasGupta (2008)), the first term is upperbounded by 0.8 (04, nj n-2 ^^0 "°Ve 
and converges to as ri2 — > 00 since inf^ogRTn Yg [£eo,ni] 7^ by assumption. To evaluate the second term, we 



just have to consider the event such that 



r„ 



< e 2 n 4 ^ since the probability that this does not occur 



converges to due to (M.S. 2) and Chebyshev's inequality. Due to Lemma[71 we can suppose that Yg [£sa nj 

-1/2 
is continuous in Bq at Oq ~ 6 without loss of generality. Therefore, Yg 



F - 



— Yg [£g,ni] 



-1/2 



< e 



for large n. Let c :~ y + n\m, aiai^m- Since 



;fe(v. 



F - 



^e[£e^n,Y-]y^i:% 



E^"«P/"" \y-ce<^Ne 



< cen 2 due to ^12\ . this is upperbounded by 






rriTI 



< $ (y + ce) - $ (y - ce) + 0.8 (a4,„i ) n^ ^Eg "" Ve 

Here the inequahty is due to Berry-Esseen bound. Letting n — )■ oo and £ — > 0, the last end converges to 
0. After all, we have our assertion. 

In TO > 2-case, the first term is evaluated using multi-dimensional version of Berry-Esseen bound (Chapter 
11 of DasGupta (2008)). The second term is evaluated by analogous but more complicated analysis. 



3.5 On logarithmic derivative and Fisher information 

Hayashi and Matsumoto (1998) gives representation of (7*5 (G^jA^) using Fisher information J^ of the clas- 
sical statistical model { P^ }«(=«■ ^'^ {Ge,M) = lim„_i.oo infivi" nTiGeJg ^^- They exploits the fact that 
the minimum variance of locally unbiased estimators equals ( J^") and achieved by T^^ g — X]"=i [^9 ^^ ) ^e^i + 
9^ , with Z^j 's denoting the logarithmic derivative. Since their regularity conditions are different from ours, we 
examine here how far this statement holds in our setting. First, we define l^^ as the Radon-Nikodym deriva- 
tive dtr5,pf"M"/dtrpf"M". Let fi^" (A) := trcr»"M" (A) (cr > 0) and p^" := d (trpf"M") /dfi^" (, 
which exists since tr cr^"M" (A) = implies Af" (A) ^ 0). Since d, / /p^"d/i'^" < (sup |/|) ||<9*/of "||p there 
is ii function d^p^" such that d^ J f p^" d^^'^ = / /9^p^"d/x'^". Using this, /^" = ^^p^" /p^" , if the RHS 
is finite. 

Lemma 8 Suppose that (M.l) and (M.2) holds. Then 1^,^ exists and sup^/|„ v^J^ vJ^ is finite. Also, if 
dipe G Tc (H) exists, and pe > 0, l^^ exists. 

Proof. Suppose (M.2) holds. Since pe > and M (•) > 0, trpf "M = means pf "M = 0. Therefore, due to 
(M.2), tid.pf'M = i (tYL'^^^pfM^trMpfU^^j = 0. Therefore, if J exists. Define Z^^ := E"i«*^^z"- 
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0< 



Let At := < w; te < 1 1^^ (w) J < (t + 1) e > , and denote by uji, the one satisfying ( Zg"^ (oj) J = te. Observe 
^ z;.t;, Y. trpf" {l^, l^^ (c.,)} A/" (AJ {i^^^^- - Z.^J K)} 

j t t 

= v^ j/'"v-25^Z,^: i^,) ^ /,^: (a.) trpf "Af« (do;) + J^ {iZ (^O) ' trpf "M" (A J 

= v^ J,^'"v - Y, [l^; K)) ' tr pf "M" ( AJ - 2^/,^: (c. J (^|^ Z,^; (c.) tr pf "M" (dc.) - Z,^; K) tr pf "M" (A,; 

< v^ Jf '"v - I [C: {CO}) ' pf-M" (d^) + e + 2E /^ l^e^'v (^0 



CM + CM 



tr pf"M" (da;) 



< v^ Jf ■' 



"v- /'(CH)%f"M"(dc.) + £ + 



2e 



,T fS^n 



which, with £ -^ 0, implies v Jg v < v J^ ' v < cxo. Also, suppose pg > 0. Then, trp^ Af = means 
AI — and tr dipf^M = 0. Therefore, the second assertion is proved. ■ 

The RHS of ^ and (JH) is lowerbounded by infM" TrGg (J^ ) , due to Schwartz's inequality. Achiev- 



ability, in fact, also holds. Define £^„ := |Mg,r^gl by 



M",L 



X{a;„;||;^J|<L}'e,j ' ye 



M"-,L 



J 



M",L 



,M",L 



9^. Obviously, f^,^ satisfies (E'). Therefore, 



^e h,i hj ' ^^°^ ^n,e ■- l^i=\ 
due to LemmalU 



IVI",L 



Therefore, E^ ^ also satisfies ([J). Also, TrG^Ve 






TrG( 



(jr-)" 



Hence, it remains to show 



lim w^J^^-^w = v^ Je^"v, Vv. This is true since v'^ Je""v = / f^" , i;»;^;) P^^ (dw„) < oo. 

Therefore, logarithmic derivative and the Fisher information can be used to represent C^ {Gg, A4). How- 
ever, it is not possible to show their chain rule, which is at the heart of the argument for one-way semi- 
classical setting in Hayashi and Matsumoto (1998). Therefore, in the next subsection, we use somewhat 
different method to prove the asymptotic Cramer-Rao type bound in the semi-classical setting. 

1 4. SEMI-CLASSICAL MEASUREMENT 

4.1 DEFINITIONS, REGULARITY CONDITIONS, AND MAIN THEOREM 

An important subclass of measurements is semi-classical measurements, which are composed adoptively in 
Rn ( < oo) rounds. At each round, we measure each sample separately, and the measurements of the r-th 
round depend on the previously obtained data. We denote by 2;,-,k(€ M') the data obtained at the r-th round 
from K-th sample, and z^ is the data (zi.i, zi.2, ■ • • , Zr^n) obtained up to the r-th round. The measurement 
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acting in the r-th round on K-th sample is denoted by Mrr^ ^. Without loss of generality, we suppose that 
in the first round the measurement is chosen deterministically. Rigorous mathematical description of such a 
process is given in the following subsection. 

Define z^,. := (z^-i, z^a, • • • Zr.w), and 'lX,k -^ (zr,K, z^^k+i, • • • ,2fl.„,n)- ^r, ^r,K, S^,, and *B|„ is the 
totality of Borel sets over the space where z^ , Zr^^, z^^, and z|^ is living in, respectively. The instrument 
corresponding to successive application of Mi,i, Mi^2, ■ ■ ■ , Mr'^^^^ and M^.'^^ \ M^''^^^, • • • , M^""^^ is denoted 
by M^ ^ and MJ,^''"\ respectively. Note that they depend on n, although we do not denote the fact explicitly 
for the sake of simplicity. Note also that i?„ is arbitrary but finite. 

Note that in other literatures such as Hayashi and Matsumoto (1999), the term 'semi-classical measure- 
ment' refers to more restricted class of measurement, which is called one-way semi-classical measurement in 
this paper. The restriction is that in r-th round, we measure r-th sample only (Hence, i?„ = n). 

The asymptotic sem^i- classical Cram,er-Rao type bound Ce {Gg,Ai) is defined by 
Ce iGe,M) := Im inf {nTrGeMSEe [£„] ; M" in "H®", semi-classical, dU, (E)} . 



Theorem 9 Suppose (M.1-3) hold. Then, 

Ce {Gg,M) = inf {TrGeYe [Soo^] ; MHn H, H, (E') with n = l} , 

= inf {TrGeYe [Soo.i] ; MHn H, H, (E) with n = l} . 

4.2 ADAPTIVE MEASUREMENT 

In this subsection, we give mathematically rigorous account on a composite measurement N M of measurement 

M followed by N", where N" is composed depending on the data w e M' from M. More specifically, 
a; ^ N"^ [A'] can be approximated by a sequence of simple functions except for w e A where M (A) = 0, so 
that the function is strongly measurable with respect to Pj^ for any p. We show this composite NM can be 
described using an instrument. (The contents of this subsection should be well-known to specialists of the 
field of measurement theory. The author, however, could not find a proper reference.) 

The key fact is Theorem4.5 of Ozawa(1985), or that there is a family of density operators {p^} p„, 
(a posteriori states) with J^^^trAp^P^ (dw) = trAM [A] (p) (VA e BiH)). Since B {H) is the dual of 
Tc{'H) with the pairing (p, A) := ii pA (TheoremII.7.2 of Holevo (1982)), Ozawa's statement is equivalent to 
the weak measurability and the existence of Pettis integral of the function w — > pj^ . As summarized in the 
end of Subsection 2.1, a; — > pj^ in fact is strongly measurable. Also, since /^^^ ll/'w lli ^T i^^) — 1 < oo, 
the Bochner integral Ji^p^P^P!f (dw) = M [A] (p) is convergent. 

First, we show Pp^ is well-defined. Since a; — > N"^ [A'] and w — > pj^ are strongly measurable, they can 
be approximated by simple functions. Therefore, a; ^^ tr N" [A'] pj^ is a measurable function for any A' e 
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NM 



<B (M'), and Pj^'^ (A x A') := /^^^ tr N" [A'] p^Pjf (duj) is well-defined and cr-additive. Therefore, P^ 
can be extended to *B (R' x M') due to Hopf's extension theorem. Moreover, with A^^ := {a;'; (w,w') G A}, 
/^ tr N'^ A(^ pj^ Pp^ ( d w) exists and equals F^^^ | AJ for any Borel set A; Let £> be the totality of A such 
that the assertion is true. Obviously, 2) is a Dynkin system, and contains cylinder sets. Therefore, due to 
Dynkin's lemma, £> = «8 (M' x R'). 



Next, we show that pj^ is well-defined. Since f 



A,, 



Pj^ (duj) < 1, the Bochner integral 



NM 



A 



(p) := J^ N"^ A„ p^ Pj^ (dui) is convergent. Also, its trace equals Pp^ {^(^ ^'^^'^'^ ^^ ^^^ / '^^^ 



be exchanged due to Fubini's theorem. 



In addition, p 



NM 



A 



(p) is affine and completely positive, as is proved in the following. Observe 



Bochner integral J N'^ 

L n-[a. 



A, 



P^ (dcj) in B{tc{H)) is well-defined in terms of ||-||^j, due to 



Piy (dcj) < 1. Therefore, there exist sequences of families -^ N^ ' > and -^ A^- > (fc 



I,---. 



oo) of completely positive maps and Borel sets, such that for any p 



NM 



A 



(p)-ENfM 



A 



(fe) 



(P) 




Aj-Y^Hf\^,.,\p^P^{du.) 



< 



A, 



y^ nC^)-. ^ -1^ 



•ij X^W > Pu, 



P;(dc.)< 



Aj-ENfx.c 



Pp^(dc.)^0, 



asfc^oo. Since ^N^^M 



A 



(k) 



is affine and completely positive, we have our assertion. 



Finally, A ^ NM [aI is an instrument. Obviously, trNM [M' x R'] {p) = 1. Also, 



NM 



Ua. 



ip) 



U A... 



P!:'^;(da;) 






A, 



P!:'^;(da;) 



E 



N" 



A 



J,w 



p::fp;(dc.)=ENM 



A, 



(p), 



where the third identity is due to Fubini's theorem of Bochner integral. Therefore, A ^ NM 
(T-additive in terms of strong operator topology in B {tc {%))■ 



A 



4.3 LEIBNIZ RULE 

In this subsection and the next, so far as no confusion is likely to arise, we denote Pg " and E^ " by Pg and 
Eg, respectively, where Mg is a semi-classical measurement. 

Lemma 10 Suppose T„ satisfies (E'). Suppose also (a) pg > 0, V0 G O, or (h) the estimator is one-way 
semi- classical. Then, 3A G ®„ s.t. M"(A) == and Eg [T„| *8^ ^] (z^ ^j *■' continuous in 9 for Vr Vk, i/ 

z„ ^ A. 

Proof. The case (b) is due to the fact that Ml^r^^ acts on pg "^ ^ . 

Eg [r„| <B^_^] (z^,,) = /T„ (z|_^,z^_i) tr pf''^'^^^^ mI^/-' (dz|_^). Here, by abuse of notation, mI^/-' (A) G 
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z^t^eA 



B ("H**" ''+^) in case of one-way semi-classical measurements. 

The case (a) is proved as follows. Let ge,6' (4.^) — ^O' [Tn\ "^i^] {4,n) - Ee [Tn\ "Bj^^] {4,n)- Observe 
U.eA^o [Tn\ si,.] (4, J Pg '^■'' (dz^,J, which equals 

/z^ eA -^" (^r-.KJ^r.fc+i) trpf"Af" ( dz^,,dz]!,,^^ j, is continuous in for any Borel set A, due to ([7]) of 
Lemma|31 Also, observe 

9e,e' (zD Pg "■" (dzi; ,,) 
zy.^eA 

/ (e,, [t„| si. J (zi. J f;'- (dzi, J - E, [r„| ^ij (zi, J f5- (dzi. J j 

Jzi^eA I J 

Ee. [T„| ^ij (zi,,) |p;'- (dzi^J - Pf- (dzi^,] 

r„ (zi:^.,z;,+i) tr (pf," - pf") Af« (dzi^^dz^.+i 

+ I E,, [T„| ®i, J (zi,,) tr (pf," - pf ") M^^, (dzi.,dz;,+i 

Tending 6*' — > 0, the last end converges to 0, and so does the left most side. Hence, due to bounded 
convergence theorem, we have j^i ^^ [limgi^g gg,gi (z^^)] Pg '''" (dz^,,) — for any Borel set A. Therefore, 
gg^g> (z^^) = for Pg-a.e. Since Pg (A) = trpeA/(A) = ^ PeM (A) ^ due to pg > 0, the proof is 
complete. ■ 

Lemma 11 Suppose (E') is satisfied. Suppose also (a) pg > 0, V0 £ 8, or (h) the estimator is one-way 
semi- classical. Then, we have Leibniz rule: 

dr^9 \Tn]\g=g, = [9,E,E,„ [r„ | ^i^ + dj^g.^g [T„| S^:, J ] ^^^^ . (15) 

P^ is a special case of ([15]). To see this, observe that the estimator (|10p is viewed as semi-classical 
considering the quantum statistical model {pf" j/jpQ- 
Proof. 



5.E, [T„]|,^,^ = a.E.Ee [r„| «i_J 1^^^^^ 



lim 



E,,E,„ [T„| »i:_J - Ee„Eeo [T„| »J:^J Ee- [E^, [r„| «i,J - Ee^ [T„| S^^J 



(16) 



where the convention is that 9'^ = 9^ {j ^ i). The first term converges to di^jg'&gg [T^IOS^^] , due 

I Ego \Tn\ ^r k\ I < ^i-n and LemmaSl Observe that the second term should converge due to the convergence 
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of (|T6|) and the first term. Moreover, 

Efl, [ Eg- [Tl I ^i^,] - E9„ [Tl I 'Bi:^,] ] Ee„ [ E^- [T^ \ S^^ J - E^^ [T^ \ ^i. 



I (E,, [T^l 33iJ - E.JT^I Sy ) tr 



9TI „Q9n 



-Sn 



<,...(d<J 



(i) 



_9_ 
9^ 



E,, [T^l »J;_J - E,„ [T^l ©i^J) tr9,pf"M,^^^,,^ (dz^^, 
< j |E,, [T^l «i:_J - E,„ [Tl\ ^i^,] I tr0lV„,„A4„,,,. (dz^^J ^^^ (0' ^ ^o) • 

Here, (i) is due to mean value theorem, where is a point between 9' and ^o- (ii) is due to 
I Eg' [T^|*B^,-] —Ego [T^|*B^,-]| < 2a4.„ and LemmalU (iii) is due to LemmalTUl Therefore, the second 
term equals 
lim ^''o[E«>[^^l'gU-E>>o[^^l'B^..]] ^ 5^Ee„Ee [n\ ^i,] e=eo- After ah, we have (HH). ■ 

9'—>-9o II '^11 ' 



4.4 ON LOGARITHMIC DERIVATIVE 

Applying Leibniz rule to the indicator function, we can prove that 



a,/^p;'-MA')d<'-(-.,K) 



is finite. However, in general, one cannot prove existence of diPg "^"^^ (A')- Therefore, we cannot define 
logarithmic derivative of the conditional probability distribution Pg ''"^^ , nor cannot use the argument in 
Hayashi and Matsumoto (1998) in semi-classical case. 

In one-way semi-classical case, which is treated in Hayashi and Matsumoto (1998), one can safely define 
the logarithmic derivative of Pg '■•''+^, since Pg "■+^■'+1 (dzr+i,r+i) — trpe^^J k+i ('^-^r+i.r+i)- Therefore, 
their argument can be made regorous, though we do not go into detail. 

4.5 PROOF OF THEOREMS 

Observe that the estimator pUj) with rii = 1 is one-way semi-classical. Therefore, the achievability by 
(one-way) semi-classical measurement follows from Lemmas[5]|6l Therefore, below we prove the lowerbound. 
In case , pe > {Vd e 6), due to the proof of lowerbound part of Theorem[Tl we have the following 
lowerbound. 

C0{G0,M) > Im inf {TiTr Ge Ye [f e,n] ; semi-classical, @, (E') } 
In the following, we reduce the optimization over semi-classical measurements to the one over independent 
semi-classical measurements, or one-way semi-classical measurements such that Ne^.^ acting on K-th sample 
cannot depend on the data y^' from Ng^^K' (^ 7^ k'). 
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Lemma 12 Suppose that semi- classical estimator £eo,n — {Teo,n, Mg } satisfies 0j, (E!'). Suppose also 
pe > 0, V0 G Q. Then, we can find an estimator J-'n.ea — {'S'so,"! ^e } ' such that N^ is independent semi- 
classical, Yog [Tea,n] ^ V^o [fffn.n], (0), and (E'J hold. Moreover, Sg^^n is in the form of |_??| ), where Fg^^,^ is 
the function of the data y^ from N^q^k, such that Egg [-F6Io,k] = ; 



Sgo,n — 2_^ P9o,K {Vk) + 6q . 



K = l 



Proof. Since Eg 



TgJ'^^ 



r,K+l 



a.Ee [Tg„^n]e^g^ - ^ \^9o 



satisfies (E'), we apply Leibniz rule (llSp recursively to obtain 

d 



E« 



d 






_d_ 

=60 ^ W 



Q0i (Ee Ego Te„^n\'Sj^^.,^_^ 



Eflo 



Eg 



Efl 



7(9o,n|S 



R„.n-1 



i'»k,«-2 



(^Ee p^o Eg„ Too^n\^R^^n_i 



l'Sk.n-2 



5:Ea^(E«o 



r=l K=l 



Eso 



■Ee 



Eg 



■Eso 



J0o,n|®fl„,ri-1 



•••|«L+1 



I Si 



■■■l«ti 



(17) 



ZlzJa^V^'^o -^'^ ^''o T'Oo,n\^r,K+l 



m 



r=l K=l 

Observe that, conditioned by S^-i, the random variable Zr^K and Z^.k' are independent, due to the 
composition of the measurement. Therefore, due to Fubini's theorem. 



Efl 



Efl 



Efl 



-teo,n|-Or,K+l 



IS 



7\K. 



ISr-l 



E0„ [re^,,,!^,] [| dPeolzr.K'ISr-OdPelz,,,!^,-!) 



= Ee [[EjIq [Tgo^nl (Sr-l, Sr.ft)] |Sr-l]] 

Therefore, 

i?„ n 



d 



S.Ee [T,„,„],^,^ = EE 9^ (E,„ [Ee [E,„ [T,,.,| (S..1, «.,«)] |»._i]] 

r— 1 K— 1 

Let us define, with the convention *8o = {0,M'}, 

feo,r,K '■= Ego [Tflo,"! (*Br-l, *Br,K)] — Eg^ [r„|*Br-l] , F0o,K '■— J2r=l fOo,r,K ■ 

Since fg^^r.K also satisfies (E'), we can apply Leibniz rule ([15]). Therefore, 



(18) 



diEgfgg^r,K\g^g„ = diEg [Eg [fgg,r,K \^r- 



lJJU=fl 



= idiEg„ [Eg [fgg,r,K l^r-l]] + diEg [Eg,, [fgo.r.K l^r-l]] )g^g„ — diEg„ [Eg [/flo.r.K |Sr-l]]le^e„ 

(i) 

= a, (E,„ [Eg [{Eg„ [Tg,,,\ (S._i, S.,.)] - Eg„ [Tg,,,\ »._i]) |S,_i] ]),^,^ 

= 5,Ee„ [Ee [Eeo [Tg,^n\ (S.-l, ».,«)] |Sr-l]]le=eo " (1^) 

(u) 

Here, (i) is due to E^,, [fg„^r.. |*8r-i] = E^,, [(E^^ [Tg„,„\ (*B.-i, «.,«)] - Ee„ [r„|<B.-i] ) |S.-i] = 0, and 
(ii) is due to Eg [Eg, [Tg.^nl'Br-i] |Sr-i] = ^Bo [Tg o,n\'Br-i]- Combining dUl) and dUl), we have a,Ee [Tg„,n]\g^g 



16 



Ef=i E:=i 9.E, [fe,,r,.]\o^Oo- Therefore, with 5^^^„ := ELi ^«o,k (^rJ + ^o, {m^„, 5^„,„} is bcahy 
biased at Oq. Also, observe the following relations: 

^eofeo,r,Kifeo,r\K') =0 {k^k! or r ^ r') , ^eo feo,r.K{Teo,n) —^Qq feo,r,K{feo,r,:i) ■ 
Due to them, the variance of this estimate is not larger than the one of Tg^^n'- 



Below, we define Ng ^ . First, using Pg , we prepare n oi fake ensembles 



® ■■■® P9o 



{k = 1,- • • ,n), composed with single pg and n — 1 of pgg. Then Ng ^ is the application of Mg to K-th fake 
ensemble. We denote by z,.''^, the data obtained at r-th round from the K'-th (possibly fake) sample in the 
K-th fake ensemble. The data y^. from Ng ^ is y^ '■— z^ . 

If 6' = 9q, Yk = 7j^' obeys the same probability distribution as Z/j^ , for any k. Therefore, Vg„ [{Ng , Fg^^^ (^k)}] 
equals Ve„ [{W^^.Fg,,, (Z/?J}]. Therefore, dueto Ve^ [{m^'^, ^^^ „}] = ELi V^o [FboA. wehave V^o [^„,eo] = 



V« 



jl^eo' '^eo,«} < Vso ['S'eo,n]- Analogously, we can also show Eg^'" [S'eo,„] = Ego 



'-'«„ 



Finally, we show SiEg"" 



Si 



- a,Ee 






= 51 . Observe 



86^ 



(e:^V.„,..(z11\,z(';)))^ 



d_ 



/.„,.,. (zi^_\,4r.on n dPg,, [zi^i^4';u) dPg {zi?y^^u] 



r' — l \k':k'^k 

= d, (Ee [Eso [• • • Eg [Ee„ [Eg [/eo,K,,r|*Br-i] | (*Br-i,K, *Br-2)] | 'Br-2] • • • l*Bi,K]])g^g„ 

= 9i (E^o [Eg [/9o,K,r|23,._i]])g^g^ + z^ ^i (Eeo [Ee [Ee,, [/eo,K,r| (*8r\K,*8r'-i)] | *Br'-2]])g^g;, 

r'=3 

+ 9i(Ee [Ego [/eo,K,r|*Bi,K,]])g^e^ 
where the last equality is due to Leibniz rule. Due to the definition of fg„.K,r, Egg [fgg^n.A {^r',K, ^r'-i) 



Si 



= ^^Eg 






— (5^ follows 



{r' <r- 1) and Eg,, [/eo,K,r|*8i^K, ] are zero. Therefore, di&g"" 

from (fT9|) . Trivially, J>i,eo = {Ng , ■S'eo.n} satisfies (E'). After all, we have the lemma. ■ 

Lemma 13 Suppose Ng is independent semi- classical. Suppose also that Sg^^n is in the form of jl'T^ , and 
that J'nfia — {Ng , Sg^^ri\ Satisfies (0j and (E'). Then, we can find an estimator £q i — {Mg , Tg -^} acting 
on single sample, with ^, (E'), andYg^ £g^ i = nYg^ [Tn.eo]- 

Proof. Mg is constructed as follows; generate x^ S {1, • • • ,n} according to uniform distribution, and apply 
N^o.^K to pg, generating the data y^. The data resulting from Mg is the pair y^ := (a;^, 2/a,„). Tg ^ is defined 

Observe Eg^"" [t^^ ^J - iELi^Eg'"" [^^o,.] + ^o = Eg^° [Sg„,r.], implying @ for 8'^^. MSE of f ^^, is 
computed as follows: V^o [f^„,i] = ^ELi"'Eg^° [f^^,, (FJ i^e„,„ (r«)^] = nVe„ [J^n,eo]- @ and (E') for 
S'a are trivial. ■ 
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Due to LemmafTOl the above two lemmas leads to '>'- part of the first identity of Theorem[9l in the case 
where pe > 0, VO G Q. The statement for the general case is straightforward consequence of the following 
lemma. 

Lemma 14 Let fe.„ := {Mg , Tg^n} be a locally unbiased estimator at 9 with (E'). Then, 3V > 0, Ve > 0, 
3£'^ g I '■— {M'g,T^ g ^} acting on a single sample with (E') and Wg [£e,ii\ > — Ve 



'-'6,1 



(1-e) 



V. 



Proof. Let cr > 0, and define pe,e := (1 ~ e) pg + ea. Denote by Eg^^ [Te.n] and Ve.e [£0,n] the aver- 
age and the variance of £g,n with respect to pe,^. Then, there is a locally unbiased estimator Eg^e.n '■= 
{Mg, Te_e,„} with respect to the quantum statistical model {pe,s}g^Q which satisfies (E) and the relation 
Tg^e,n — Y^ {Te,n ~ ^B,e [Tg^n]) +&■ Sincc pg^^ > and the family {p0,e}g^Q satisfies (M.1,2), Lemmas lT^lfOl 
imply existence of a locally unbiased estimator £'^ g i '■— {Mg,T^ g ^} acting on single sample such that (E') 
and Ve,£ [£e,e,n] > ■k^e,e £g.e,i ■ 

Defining £^1 := {Mg, T^,J by T^ i := (1 - e) fe,e,i -Eg T'g^^ ^6*, £:^_i is locally unbiased with respect 
to {pg}g^Q. Letting V be the variance of £g_n with respect to a, we have 

Vo [£e,n] = (1 - e)Yg,e [£e,e,n] " y^V^ > ^^V0,e [^^l] ~ T^^ ^ ^^^ [£'g^^ - ^V . M 

Combining with the achievability, the first identity of Theorem 12] is proved. The second identity is shown 
using the analogous argument as the proof of (l3l) = (|4]), and Theorem[9]is proved. 

Note that the estimator achieving the lowerbound is one-way semi-classical. This means that the optimal 
asymptotic cost of one-way semi-classical measurements is also Cg {Gg,A4). 

5. LOCC STATE ESTIMATION 

Suppose the state pf"" is shared by remote party, Alice and Bob: V, = Tia ® T~ib- Suppose also that Alice 
and Bob can exchange classical messages, but cannot interact quantumly with each other. So they do the 
following i?„-round measurement M": At each round, Alice and Bob measures her/his share of the samples, 
and the measurements of the r-th round depend on the previously obtained data. We denote by ^"(e M') 
and ^^(e M') the data obtained at the r-th round by Alice's and Bob's measurement, respectively. Also, ^^ 
denotes (Cr';'?r')r'=i- ^^ denote by C^ and £J: {x — a, b) the Borel field over the space which ^^-i and ^^ 
takes values in, respectively. The measurement acting in the r-th round by Alice (Bob) is denoted by M^^j^^ 
( M^''^~\ resp.). Such operations are said to be local operations and quantum communications {LOCC, in 
short). An important subclass of LOCC is local operation (LO), where Alice and Bob does M"^^ and MJJ, 
independently. 

The difference between semi-classical and LOCC measurements is the split of the actions. In the former, 
split is between samples. In the latter, the split is between Alice and Bob. Other than this point, basically 
they are the same concept. Especially, LO corresponds to independent semi-classical measurements. 
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Define Cg ' {Ge,M) and Cg {Gg, Ai) by restricting the range of measurement to LOCC in the definition 
of Cg {Ge,M) and Cg {Gg,M), respectively. Then, trivially, we have: 

Theorem 15 Suppose (M.1-3) hold. Then, 

Gl^ {Gg,M) =ini{i:TGg„\g,[£g,s] \ LOCC, m, (E) fornix}, 

Gf''{Cg,M)= lim inf{nTrGeoV(,„[£eo,n] ; LOCC, 0), (E)}. 

?i— >-oo 

From here, we focus on the case of pg ~ P% ® Pg- The motivation of studying this seemingly easy 
case is as follows. Suppose that rankpe ~ 1, diniH < oo, and pg ^ p^ ® p^. Then it is known that 
Cg — Cg ~ Cg — Cg (Mat suffioto (200 7) ) . The estimator used to show the identity, however, fails in 
case of pg = Pg (^ Pg. 

In this case, we can translate the argument in Section 4 about the lowerbound to the asymptotic cost of 
semi-classical estimators to LOCC estimators. To see this, observe the proof of Lemmas [TTI - [T2l are valid even 
if non-identical independent samples {2)"=i Pg^ ^i'g given. Below, we present the analogue of Lemma [T2] 
(The proof is omitted being almost parallel.) Using this lemma, the optimization over LOCC is reduced to 
the one over LO, where Ng ^ and Ng ^ is are measured independently, producing the data y° ,y^ . 

Lemma 16 Suppose pg = pg (^ pg > 0, \f9 G Q. Suppose an LOCC estimator Sg^^i — {Mg , Tg^^n} satisfies 
(0) and (E'). Then, we can find an LO estimator Tg„,n — {Nep.n, Sg^^n} , such thatYg^ [Tgo^n] < V^q [£go,n], 
(0), and (E') hold. Moreover, Sg^^n 'is in the following form: 

So,,n - Fl^,, {^D + Fl^,, (eO + 00 . (20) 

Let M.X denote {pglggg, [x — a,b). Observe that the map 9 -^ pg may not be injective. Therefore, the 

vector space {vg ; X^I^li v^'^diPg = O} may not be {0}. We denote by Hg the projector onto the orthogonal 

complement of this vector space in M™. Letting J"^^^ „ := |N^^ „, F^^J^, Bg„ J"g^^ „ v^^ = if n^„v^^ = 0. 

Therefore, there is a matrix Wg with Bg^ Tg „ = Wg Hg . 

We want to minimize the variance of locally unbiased estimator in the form of (|20p . First, for a given 
(^e"o.n, <,„)> we define Sg,,^ [^e"o.n,<,„] := K,,r.Fl,Ail)+ ^B,,nFLn {^n) +^0, whcrc (a,"^^„, <_„) 
moves over all the mxm real invertible matrices. Elementary but tedious calculation shows that the variance 
of such estimators is larger (in the sense that the difference is positive semi-definite) than 

K (v,o Kjy'Ko+ni {Ye„ Kjy'iiiy\ (21) 

where S^^^ := {n|^ „, ^^^^ „}, F^^„ := (W^e^J~Ve'o,n and {-y^ in ^ denotes generalized inverse. 
Observe that £g n {x = a,b) satisfies 

Bgy£l^^]=ni (x = a,b). (22) 
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(CM.l) dike, didjAe exits and are locally uniformly continuous, ll^jAeH^f, , \\didjAg\\^^ < ^ < oo, V6' G 9 



(CE) / ||r„ (w) - e\\ t^w^ (D,,, , E„) M" (dw„) < naia4,„, / ||r„ (w) - 9^ ti P^- (dw„) pf " < n (a4,„)", V0 G 6 



Table 2: Regularity conditions on models and estimators in operation estimation. 

In the end, we optimize (PT|) with the constrain ^ and (E'). Here, constrain (E') can be replaced by 
(E) without increasing the infimum, due to the analogous argument as the proof of ([3| = (|4l). 

So far we had assumed pg > Q. However, this assumption can be removed due to an analogue of Lemma fMl 
(the proof is omitted, being straightforward.). Therefore, letting Vg „ {M.) the totality of Ve [£e,n] with (f22|) 
and (E) ( x = a, b), we have 

Theorem 17 Suppose pe = Pg <^ Pg, '^0 e O. Then, 

C^ {Gg,M) > inf JTrGe (HgK"'ng + H^y-^Hg)"' , V, G VI, (M) , x ^ a, b} , 

C^-"- {Gg,M) > lim inf \nTrGg {U^gV-^U"g + H^V^-^H^)"' , F, G V,\ (M) , x ^ a, b] . 
The achievability of the lowerbound is also true. This theorem leads to a necessary and sufficient condition 
for C^ {Gg,M) {Cg {Gg,M) ) to equal C^'^ {Gg,M) {C^ iGg,M), resp.). These topics, however, wih be 
discussed elsewhere. 

6. ESTIMATION OF QUANTUM OPERATIONS 

Suppose we are given a family of completely positive and trace preserving maps C := {Ag}g^Q. Here, 
Ag : B (H) -^ B{H'), 9 ^ Q, and Q is an open region in R™. Our purpose is to estimate 9, by measuring 
the output of Ag after sending the input state for n times through it. 

Our input state p" is living in 'H^/C, where dim/C is arbitrarily large. /C may be used to store the input 
state before and/or after application of Ag . Between the K-th and {k + l)-th use of Ag, one can apply an 
operation S" : B (H' ® JC) ^ B [Ti. ® ]C) . S" may be a composition of measurement followed by preparation 
of the state to be send through Ag. After n times of use of Ag, we obtain n"=i { (^e ® I) o S" } (p"). We 
measure this by M", obtaining the data a;„ G M'", and compute the estimate T„ (aj„). The pair £„ := 
{p", {S"}"^j^ , M",r„} (, or sometimes sequence {£n}^^i also, ) is called an estimator. The probability 
distribution of the data is Pg^ {w„ G A} = tr Af" (A) HLo K^e ® I) ° -"} (p")- 

Regularity conditions, other than ([T]) on estimators, are listed in Tabled Note that they are honest 
analogue of (M.l) and (E). In the table, convergence is always in terms of IHl^^, and Di^g : B {'H) — > B {V.) is 
an affine map with D^.g (E)I{p) > diAg^ ® I [p) (||0o — 9\\ < 02), whose existence is certified by (CM.l) and 
an analogue of Lemma[3l Also, ro" (D^^e , E„) is defined by replacing diAg in di n"=o { i^s ^ I) ° S" } (/?") 
by Di.e. 
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Define C^'°^(G'0,/:) by replacing (E) in the definition of C^'°P(Ge,X) by (CE). Tlien we liave, 
honestly modifying the argument in Section 3 (the proof is omitted), 

Theorem 18 If (CM.l) holds, Cf^^iGe^C) > lim inf {nTr Gg Vg [f e «] ; ®, (E')} 

Also, the achievability of the lowerbound can be proved, with some additional regularity conditions. It 
is known that in case of Ag (/?) = UopUg, with Ugllg ~ UgUg = 1, there is an asymptotically unbiased 
estimator with TiGgYg [£„] — O (^) . Therefore, Cg {Gg,C) = for such models. In this subsection, we 
show that such a phenomena can occur only at the surface of the space of the quantum operations. 

Theorem 19 Suppose (CM. 1-2) and (CE) holds. Moreover, suppose Be > s.t., h.Q + 'Y^^^iu'^diKg is 
completely positive, for^u : \\u\\ < e, V0 G 6. Then, Gg {Ge,C) ^ 0. 



Proof. Define Aa 



Afl 



Eti ^ (^* - ^o), ^o„ ■■= {Ae„e; Et 1 1^* - ^ol < £}• Suppose Eg, 



{p^^^, js^Jg 1 , M^^_^,Tgg^n} satisfies (CE) and © at 6* = 6*0 as an estimator of C := {Ag} 



eee 



Tg„,rA'^n)tT 



dkg 



Ml (d a;„) n A,„ + ^ ^ (0' - 0^) ® I o S;' 



t=0 



,eo 



e=eo 



Here, the first identity holds since (CE) implies an analogue of Lemma|4l Therefore, since l^gafla = Ag,,, 
Edo,n satisfies ^ sX Q = Qq s& an estimator of Cg^, also. Moreover, the variance of Eg^^n as an estimator of 
>Ce„ and C coincide at = ^o- Therefore, it suffices to show the statement for the quantum operation model 
£e„. From here, we follows the same line of argument as Hayashi (2003) and Zhengfeng Ji, et. al. (2006). 

Let Qr^ be the a;-th extreme point of the convex region X)i=i 1^* ~ ^ol — S- Then, there is pg^^.g {x) such 
that kg^fi = J2l=iPeo,e (x) Agg^g^, J2l=iPeo,e i^) = 1, and pg„M{x) is linear in 9. Consider the family 
of multinomial probability distributions { pgg_g {■)}g^Q. The key observation is that Ag^-^^g is equivalent to 
random application of Ag„_g^, where x is sampled according to peo,e (■)■ Following this observation, given 
a locally unbiased estimator Sg^^n of Cg^ at 9 — 9o, one can compose a locally unbiased estimator of the 
statistical model {peo,e (Olgge ■ Prepare a quantum state pg , apply a sequence of quantum operations 
n"=i 1 {-^Sofi^o^ ® ■"■) ° '^K So f ' ■^here Xk. ~ P0o,e {i^ = !>■ • • ^n), measure the output state by Mg , and compute 
Tea,n- Therefore, due to classical estimation theory, the variance of Ega^n is O (!/"■)■ ■ 
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